Kohn-Luttinger pseudo-pairing in a two-dimensional Fermi-liquid 
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We consider possible superconducting instabilities in a two-dimensional Fermi system with short- 
ranged repulsive interactions between electrons. The possibility of an unusual superconducting 
paring due to the Kohn-Luttinger mechanism is examined. The quasiparticle scattering amplitude 
is shown to possess an attractive harmonic in second-order perturbation theory for finite values of 
the energy transfer. The corresponding singularity in the pairing vertex leads to a superconducting 
pairing of the electron excitations with finite energies. We identify the energy transfer in the Cooper 
channel as the binding energy of the excited pair. At low enough temperatures, the Fermi system 
is a mixture of normal electron excitations and fluctuating d-wave Cooper pairs possessing a finite 
gap. 
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I. INTRODUCTION 

Superconductivity induced by mechanisms other than 
electron-phonon interactions has been of long-standing 
interest. Throughout the last decade there has been con- 
tinuing theoretical search for unconventional supercon- 
ductivity mechanisms, particularly in two-dimensional 
systems. This interest has been, indeed, motivated 
by novel superconducting materials such as high-T c 
cuprates, organic superconductors as well as by the stud- 
ies of 3 Hc films. Currently there is no full understanding 
of the physical processes responsible for the pairing in 
those systems. 

Kohn-Luttinger effect is one of the oldest as well as 
among the most appealing and elegant physical effects, 
which might be considered within this quest. Back 
in 1965, Kohn and Luttinger 1 showed that any three- 
dimensional electron system with repulsive interactions 
between particles was unstable against a superconduct- 
ing transition at extremely low temperatures. The origin 
of the effect is that the screening of the bare interac- 
tion leads to the well-known Friedel oscillations in the 
electron density and to similar oscillations in the scatter- 
ing amplitude. The renormalized interaction acquires a 
long-ranged oscillatory component. Thus, there appear 
some regions where the effective interaction is attractive. 
This leads to the formation of Cooper pairs with non- 
zero orbital momenta / ^ 0. However, straightforward 
calculations showed that the transition temperature was 
extremely low (the estimate of Kohn and Luttingesi was 
T c ~ 10~ 40 K for some realistic parameters of the fermion 
system) . This extreme low value of T c was one of the rea- 
sons why the effect has not been much studied in recent 
years. 

In the early nineties, Kagan and collaborators ob- 
tained a number of interesting results within the Kohn- 
Luttinger theory 2 (such as cascade transitions, Kohn- 
Luttinger effect in a three dimensional system with long- 
ranged Coulomb interaction, Kohn-Luttinger supercon- 
ductivity in the Hubbard model, etc.). One of the in- 
teresting results was that the temperature of the super- 



conducting transition derived in the pioneering paper- 
was shown to be underestimated due to the unjustified 
extrapolation in the expression valid for large orbital mo- 
menta down to the value 1 = 1. The transition temper- 
atures calculated in Ref. were higher than the original 
estimate but still too low to attract much attention. 

One of the most natural issues to be explored has been 
the status of the Kohn-Luttinger theory in two dimen- 
sions. First of all, the Kohn-Luttinger physics is about 
the formation of bound states. It is very natural to expect 
that in the lower dimensionality it is easier to form bound 
states {i.e., Cooper pairs). However, a simple calculation 
of the polarization operator leads to the disappointing 
result: no singularity exists in second-order perturbation 
theory. Namely, the polarization operator reads (we use 
units h = c = 1 throughout the paper) : 



n(q) 



s[l - (2k F /q) 2 



if q < 2k F ; 
if q > 2k F , 



(1) 



where v = m/ (2tt) is the density of states at the Fermi- 
line, k F is the Fermi momentum, and q is the momentum 
transfer in the Cooper channel (q = 2pp sin <f>/2, where <f> 
is the scattering angle). Let us remember that the attrac- 
tive harmonics in the scattering amplitude in the three- 
dimensional case comes from the well-known logarithmic 
Kohn's singularity n s i ng (</>) = (1 + cos</>) In (1 + cos <j>) 
which exists in 3D on both sides of the Fermi-surface. As 
can be seen from Eq.JQ), the singularity in two dimen- 
sions is one-sided which suggests that no straightforward 
Kohn-Luttinger effect should exist in 2D. 

In 1993, Chubukovi showed that this simple scenario 
was not the complete story in two dimensions. A two- 
sided singularity exists, but to find it one should go be- 
yond second-order perturbation theory. The correspond- 
ing transition temperature derived by Chubukov reads: 
T c (l) oc exp [— / 2 /2/q] , where /o is the dimensionless s- 
wave scattering amplitude. Having applied this result to 
a realistic experiment on 3 He — 4 He mixture films, the 
numerical value was found as T c (l = 1) = 10~ 4 K. 
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Let us also mention a recent paper of Guinea et al£ 
in which the Kohn-Luttinger physics was phenomenologi- 
cally incorporated in a model of high-T c cuprates. Within 
this model the shape of the gap anisotropy has been ex- 
plored as a function of doping. 

The main idea of the present paper is to search for 
an effective attractive interaction by taking into account 
the frequency dependence of the polarization operator, is- 
tead of going into higher order perturbation theory. The 
account for dynamical screening, as we shall see below, 
yields a two-sided singularity. Thus, we are looking for a 
dynamical Kohn-Luttinger effect rather that the original 
static pairing problem as in Refs. and 0. Due to 
the energy dependence of the effective electron-electron 
coupling, the Cooper problem turns into an integral equa- 
tion, similar to the Eliashberg equation in the strong cou- 
pling theory of superconductivity^ 

Our paper is structured as follows: 
In Sec. II, we rederive the expression for the polariza- 
tion operator as a function of momentum q and Mat- 
subara frequency u). Using this result, we formulate 
the Cooper problem and derive the corresponding Bethe- 
Salpeter equation for the pairing vertex T(q;e,e'). 

In Sec. Ill, we consider spherical harmonics of the 
effective interaction Vi(lu) and show that d-harmonic, 
which corresponds to the orbital momentum 1 = 2, yields 
the strongest effective attraction. 

In Sec. IV, we use the explicit expression for the inter- 
action in the d channel and derive an integral equation for 
the pairing vertex. Studying this equation, we show that 
the pairing vertex may diverge if the incoming particles 
have high enough energies. We estimate the tempera- 
ture at which the pairing with the typical binding energy 
of ui commences. We conclude that at low enough tem- 
perature the system is a mixture of low-lying electron 
excitations and fluctuating Cooper pairs. We estimate 
the temperature T» at which the effect of this fluctuaing 
pairs becomes essential and may strongly change trans- 
port and thermodynamic properties of the system. 

In Sec. V, we briefly discuss the case of long-ranged 
Coulomb interactions. We argue that Kohn-Luttinger 
physics strongly depends on the screening properties. In 
a purely two-dimensional system we do not expect any 
superconducting instability to survive. If transport is 
two-dimensional but screening is three-dimensional, the 
system is qualitatively described by our theory and Ref. 



II. COOPER PROBLEM 

Let us start with calculating the effective electron- 
electron interaction V(q, uj) (where uj = 2nnT is the 
bosonic Matsubara frequency). In second-order pertur- 
bation theory there are four diagrams to be considered, 
which are shown in Fig. 1. If the bare potential A(q) is 
short-ranged the diagrams "b" , "c" , and "d" cancel each 
other out and only "e" contributes to the renormalized 



interaction. The latter diagram is functionally identi- 
cal to "b" but depends on p + p' rather then on p - p', 
where p = (p,s). Thus, knowing the two-dimensional 
polarization operator we readily obtain the total effec- 
tive electron-electron coupling. 

The polarization operator is defined as 



e n - 



(2) 



where Q e (k) = [ie — £k) is the Matsubara Green func- 
tion, £ k = (k 2 — kp) /2m, and e n = (2n + \)nT is the 
Fermionic Matsubara frequency. After the sum over e n 
is evaluated, Eq.(J2| takes on the form 



7T(q,w m ) = 2 Re 



/(k) 



d 2 k 



(2tt) 2 e(k) -e(k-q) 



(3) 

where /(k) is the Fermi distribution. At not very high 
temperatures T <C £f, it can be written as /(k) = 
6 (kp — |k|) and after a straightforward calculation we 
obtain: 



7T(z) 



'Re 



1 



1 



Rez 



1 



(4) 



where we have introduced the complex variable z for com- 
pactness: 



i \w ri 



2fa 



v F q 



and Vf = k^/m, is the Fermi- velocity. One can easily 
check that Eq.® reproduces Eq.Q if uj m = 0. To get 
the expression for the polarization operator II(q, uj) as 
a function of the real frequency^ one has to do the an- 
alytical continuation in Eq.Q). Let us note here that 
only the real part of the polarization operator renormal- 
izes the scattering amplitude. The imaginary part is not 
relevant to this renormalization. The latter quantity is 
proportional to the density of the electron- hole pairs. 

Let us now formulate the Cooper problem for the case 
under consideration. We are looking for a singularity in 
the Cooper channel (see the diagrammatic equation in 
Fig. 2). After averaging over the spin indices, the Bcthc- 
Salpeter equation can be written as: 

T(q;e,e') = V(q,^)-T£ T (k - p; £ , C) (5) 

xg c (k)g_ c (-k)v(k- P ',c-0, 

where (, e, and e' are fermionic Matsubara frequencies, 
q = p — p', and u> = e — e'. Let us emphasize that 
V(q, ui) is the renormalized interaction which depends on 
momentum and energy transfer. 
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FIG. 1: Renormalization of the scattering amplitude by Friedel oscillations. Kohn-Luttinger theory. If the bare coupling is 
q-independent diagrams "b", "c", and "d" cancel each other out. 
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FIG. 2: The Bethe-Salpeter equation for the irreducible vertex F (Cooper problem). 



Let us now consider only electrons in the very vicinity 
of the Fermi surface so that q = 2fcp sin (</>/2), where </> is 
the scattering angle. Following the standard route in the 
Kohn-Luttinger theory, we expand V and T in series of 
the normalized eigenfunctions of the angular momentum: 



V(q 3 w)=2^(w)*iW 



(6) 



analytical continuation in Eq.©. In the limit T — » 0, 
the procedure of the analytical continuation reduces to 
the simple Feynman rotation and all the Matsubara sums 
involved may be replaced by the corresponding integrals 
with the temperature serving as a "low-energy cut-off." 
The main result we are deriving in the present paper can 
be noticed in this limit as well. 



and 



(7) 



III. EFFECTIVE ATTRACTION IN THE 
d-CHANNEL 



where 



u<j> 



Then, Eq.Q takes on the form: 

% (e, e') = V«M - T£)7i (s, () C(C) V ; (C - e% (8) 
C 

where as usual C(C) is the Cooperon, which is the source 
of the BCS logarithm: 



(9) 



Let us note that Eq.© is exact at any temperature. 
However, we shall consider only the case of low temper- 
atures to avoid technical difficulties connected with the 



The next step is to evaluate the spherical harmonics 
of the renormalized interaction. At this point, let us as- 
sume that the initial electron-electron interaction is de- 
fined only for the energies smaller than some threshold 
value (D <C £f, which will be serving as the high-energy 
cut-off (just as the Debye frequency in the classical weak- 
coupling BCS theory). In this case when performing ac- 
tual calculations we can expand on lo/sf- 

The Z-harmonics of the polarization operator JIJ can 
be written as: 



2tt 



IT (<j), uj m ) cos 1(f) d(f>. 



(10) 



Keeping in mind that lo <C £f and evaluating the integral 
with the logarithmic accuracy, we obtain for even orbital 
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momenta I = 2n: 



and 



2 M I 3 2e F 



(11) 



3tt 



^ ( n + - 



and for I = 2n + 1 



/2 M I 1 , 2e F 



(12) 



-2n [C + V (n + 1)] 

Z£ F 



where ^ is the logarithmic derivative of the Gamma- 
function and C sa 0.577 is the Euler's constant. 

From Eas. (|llll2fl we see that the dependence on the 
orbital momentum is very weak. The effective interaction 
can be written as 



Vj(w) = TTi{lu) {a(0) (-1)"' + 2 [A(0)A(2fc F ) - A 2 (2fc F )] } 

(13) 

where A(q) is the Fourier-component of the bare interac- 
tion potential. 

If the initial interaction is q- independent, we see that 
the effective interaction is attractive only for the even 
values of the orbital momentum I = 2n ^ 0. The effective 
attraction is the strongest for 1 = 2. The corresponding 
eHiarmonics reads £ 



Vd(w) 



3 2 M 2e F 
= v\ - — In- — -. 

2tt 2e F U 



(14) 



IV. PAIRING AT FINITE ENERGIES 



We can substitute result (|14f> into the Bethe-Salpeter 
equation JSJ which turns into an integral equation (at 
T — > 0) with a well defined kernel K(e,e') = Vd(e — 
e')C(e'). One can easily see that if the incoming particles 
have zero energies, the Cooper singularity gets canceled. 
However, at finite energies the Cooper logarithm survives 
being cut-off by the energy transfer. 

For further treatment, let us define the following aux- 
iliary dimcnsionlcss variables and functions: 

x = e/2e F , 
x = uj/ef, 



(27T) 



3/2 



In these notations, Eq.© takes on the form 

g(x,x') = -g (x- x') + k [ dy 90 X V g(y,x'). (15) 

J \y\ 

The integral in (|15[) is defined in such a way that the 
large-y singularities are cut-off by cD/sp and low-y singu- 
larities at r = T/2e F . 

It is hard to solve Eq. l|15H exactly. However, we are 
mostly interested not in the detailed solution but in the 
possibility of a singularity in the pairing vertex g(x, x') 
which would be a signal of a superconducting pairing 
(but not necessarily a global superconducting instabil- 
ity). Let us emphasize here that g(0, 0) = by the con- 
struction and it can not diverge simply because there is 
no attraction in this case, unless we take into account 
the higher order diagrams. At finite energy transfers, 
the large Cooper logarithm appears which yields a diver- 
gence of g(x, x') which we interpret as an appearance of 
fluctuating Cooper pairs built up of the electronic exci- 
tations with finite energies. One of the ways to search 
for the singularity is to consider the eigenvalue problem 
for the kernel of integral equation (|15|) : 



A(x) = K j 



\x-y\ 



111 



\x - y\ 



A(y)dy. (16) 



The singularity exists if there is a non-trivial solution of 
this equation. To get some qualitative estimates let us 
approximate the corresponding eigenvector by the follow- 
ing trial function: 



A(x) = A + Ai \x\ 



(17) 



where Ao and Ai are some weak (logarithmic) functions 
of x. From Eas. i|ltjll7[l we can derive the self-consistency 
equation which yields the estimate for the threshold tem- 
perature at which the pairing with the typical energy 
transfer of u> commences i£ 



T p (u>) ~ u> exp 



1 



k 2 x 2 In (2e F /uj) 



(18) 



This estimate can be alternatively derived by considering 
the resolvent of the integral equation straightforwardly. 
Namely, one can formally re- write Ea. (|15fl as follows: 

g = go + KKg, 



g (x) = - Nl m T3' 



g(x,x') = 



2ir 



T(e,e') 



where K is the operator with the kernel in x-representa- 
tion being equal to K (x, y) 



The solu- 



\v\ 

tion of this equation has can be formally written as 



g = R(k) g a = 1 - kK 



go, 



5 




FIG. 3: Aslamazov-Larkin contribution to the conductivity. 
Small wavy lines correspond to the factor ev. Shaded boxes 
are the pairing vertexes which in the case under consideration 
are functions of the four variables F(fci, k[, k' 2 ). 



where R(n) is the resolvent, which can be also written 
as: 



R(k) = J2 K n K r ' 



(19) 



71=0 



where K n can be found by evaluating the convolution of 
the corresponding kernels in the ^-representation: 

KW(x,v)= I K(x,z)K ( - n - 1 \z,y)dz. 



Studying the geometric series (|19|l . one can see that its 
2n's term contains the logarithm ln ra (ui/T), with u> being 
the typical energy of the electrons in the Cooper channel. 
Summing up the series, we reproduce Eq. l|18|) . 

The integral equation (|15|l and the corresponding 
eigenproblem Ijlbfl are mathematically well-defined for 
an Y x T (*- e - w J$ T). However, it does not make 
too much sense to study the structure of the solutions at 
such energies in the framework of our formalism based 
on the Matsubara technique. Thus, result l|18l) has the 
following domain of applicability: 

T < uj < uX e F . 

Working in this domain, the replacement of the Matsub- 
ara sums by the integrals is legitimate and our interpre- 
tation of w > T as a real energy of a pair is valid as 
well. 

Let us now briefly discuss how the appearance of the 
fluctuating pairs affects the physical properties of the sys- 
tem. The correction to the conductivity is described by 
the diagrams similar to the ones in the conventional fluc- 
tuation theory^ (see e.g. Fig. 3 where the Aslamazov- 
Larkin- like diagram is shown). It is a rather difficult 
problem to calculate the corresponding contributions in 
the case under consideration. However, we can get some 
qualitative insight by noting that the analytical contin- 
uation to the real frequencies in the expression for the 
conductivity contains the factor coth(^), which is ba- 
sically the Bose-distribution for the fluctuating Cooper 
pairs (the density of the Cooper pairs). This factor and 
the corresponding correction are exponentially small un- 
less there exist Cooper pairs with u> ~ T. Using Ea. (ll8|l . 
we can estimate the temperature T» at which such pairs 



appear. It is defined by the condition T p (T* 
we readily obtain t£ 



e f exp 



(27T) 



3/2 



3tt 



£F 



(Ai/) 



T*. Thus, 



(20) 



At this temperature, contribution to the conductivity due 
to the preformed Cooper pairs may become comparable 
to the Drude conductivity of a normal metal. 



V. LONG-RANGE COULOMB INTERACTION 

Until now, we have been studying a Fermi-system with 
the bare electron-electron coupling being short-ranged. 
It is worth considering the case when the initial interac- 
tion is the long-range Coulomb repulsion. In this case 
our treatment is not applicable since the momentum- 
dependence of the Coulomb interaction becomes crucial. 
However, we can get some qualitative insight into the 
problem without cumbersome calculations. There are 
several possibilities one can consider: 

First, we can study a system in which both transport 
and screening are two-dimensional. In this case, we can 
readily conclude that there is no possibility for Kohn- 
Luttinger pairing because the long-wavelength Thomas- 
Fermi screening is weak 



V(r) 



d 2 q 2ire 2 1 



(2tt) 2 q e(q) 



' qr cx -4 



(21) 



where ao = 1/me 2 is the effective 2D screening length 
and the Thomas-Fermi dielectric function has the stan- 
dard long- wavelength form: 



e(g) = l + — . 

a q 



(22) 



We can now calculate the spherical harmonics of the 
screened Coulomb interaction 



Vi 



2 f 2ne 2 A 

; — cos IS dA, g = 2fcp sin — , (23) 

g + 2/ao " 2 



which are certainly all repulsive and remain repulsive 
even after Friedel oscillations are taken into account (see 
Fig.l). Thus, even going beyond this long- wavelength 
Thomas-Fermi analysis we do not expect any pairing in- 
stability for the long-ranged Coulomb interaction to ap- 
pear, as long as both transport and screening are two- 
dimensional. However, the account for the dynamically 
screened Coulomb interaction may lead to other impor- 
tant effects such as renormalization of the Fermi-liquid 
parameters (effective mass, (/-factor, etc.). This issue is 
currently being investigated by the authors and the re- 
sults will be reported elsewhere. 

Second, one can consider a system in which transport 
is two-dimensional but screening is three-dimensional. In 
this case the Coulomb interaction is well screened and de- 
cays exponentially at large distances V(r) cx exp (— r/d) 
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(d is the screening length). In the limit kpd <C 1, the po- 
tential becomes effectively short ranged and, thus, the 
theory developed in the present paper is qualitatively 
valid. Let us note that in this model the high-energy 
cut-off is basically the Fermi energy which violates the 
assumption Co <C e we used in our calculations. This, 
however, should not change main qualitative result of 
the paper. 

There is also an intermediate situation which may ex- 
ist when the two-dimensional Fcrmi-liquid lives in the 
very close vicinity of a metallic substrate. In this sit- 
uation, each two-dimensional electron produces an im- 
age in the metallic substrate so that the bare electron- 
electron interaction decays only as r~ 3 at large distances. 
In this case, there is no simple answer whether the Kohn- 
Luttinger pairing exists or not. Presumably, the Kohn- 
Luttinger pairing in such a setup is possible if the Fermi- 
liquid is dilute enough, so that Friedel oscillations may 
compete with the initial dipole-dipole coupling. 

VI. CONCLUSION 

Before concluding, we point out that earlier theoretical 
work in the literature has considered^ the possibility of 
bound states and Cooper pairing in a dilute 2D system of 
fcrmions interacting via a short-ranged repulsive interac- 
tion. Engelbrccht and Randeria have considered a regular 
expansion in the T-matrix in two-dimensions analogous 
to the expansion on the dilute gas parameter k-pa <C 1 
in 3D>ii Apart from the three-dimensional result, they 
have found an unusual pole in the particle-particle chan- 



nel. Although we do not find any obvious connection be- 
tween our microscopic analysis and this earlier work^* 
the claim of a new 2D collective mode interpreted as a 
bound excitation of two holes is somewhat reminiscent 
of our finding in this paper that Kohn-Luttinger type 
superconducting pairing is possible at finite excitation 
energies. Whether there is a deep connection between 
our work and the earlier results^ remains unclear at this 
stage. 

Summarizing, we have shown that a clean two-dimen- 
sional Fermi-system with a short-range repulsive interac- 
tion between electrons becomes unstable against a forma- 
tion of d-wave Cooper pairs with a finite binding energy 
at a low enough temperature. Thus, the low-temperature 
state of the system is a mixture of low-lying electron ex- 
citations and preformed fluctuating Cooper pairs. The 
new type of carriers may noticeably change the physical 
properties of the system such as conductivity, suscepti- 
bility, etc. at a temperature T* [see Ea.()2(Jfl]. Let us 
note that from our theory it follows that the fluctuat- 
ing pairs appear within the normal state having a finite 
gap which is connected with the binding energy. Note 
that there is no global superconductivity specifically pre- 
dicted in our theory, only a pseudo-pairing at finite exci- 
tation energies. The results obtained in the present paper 
may be relevant to the pseudogap experiments in high-T c 
superconductors^ 
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manuscript. 
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